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Abstract 

Using a five dimensional AdS soliton in an Einstein- Yang-Mills theory with SU(2) 
gauge group we study p-wave holographic insulator/superconductor phase transition. 
To explore the phase structure of the model we consider the system in the probe limit 
as well as fully back reacted solutions. We will also study zero temperature limit of 
the p-wave holographic superconductor in four dimensions. 
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1 Introduction 



The application of AdS/CFT correspondence [T] in condensed matter physics has been at- 
tracted lots of attentions in recent years. In particular it has provided gravitational descrip- 
tions for systems which exhibit superconductor/superfluid (2|]4) phases. Since in condensed 
matter physics we are typically dealing with systems at finite charge and temperature, in 
the context of the AdS/CFT correspondence the dual gravity descriptions should be given 
by gravitational models which admit charged black holes as vacuum solutions. 

Indeed the simplest model may be provided by an Einstein-Maxwell theory coupled to a 
charged scalar field. For this model, it has been shown that the charged black holes become 
unstable to develop scalar hair for sufficiently low temperature |2] and moreover the U(l) 
symmetry is broken near the black holes horizon. Within the framework of the AdS/CFT 
correspondence, the charged scalar field corresponds to an operator which carries the charge 
of the global £7(1) symmetry. Having non-zero hair corresponds to the fact that the dual 
operator has non-zero expectation value and the global U(l) symmetry is broken as well. 
This phenomena may be interpreted as a second order phase transition between conductor 
and superconductor phases. This interpretation has been supported by making use of the 
behavior of AC conductivity in these phases. 

Using a five dimensional AdS soliton in an Einstein-Maxwell-charged scalar field the 
authors of [5] have constructed a model describing an insulator/superconductor phase tran- 
sition at zero temperature. Actually since in this model the normal phase is described by an 
AdS soliton where the system exhibits mass gap |6j the dual field theory is in an insulator 
phase. On the other hand for sufficiently large chemical potential the AdS soliton becomes 
unstable to forming scalar hair which corresponds to the fact that the dual theory is in a 
superconductor phase. Holographic insulator /superconductor system with back reactions 
has also been studied in [7]. 

The aim of the present article is to extend the holographic insulator/superconductor for 
the case where the AdS soliton develops a vector hair. In fact holographic superconductors 
with vector hair, known as p-wave holographic superconductors, have been first studied in [8] 
and further explored in (9,10 . The simplest example of p-wave holographic superconductors 



may be provided by an Einstein- Yang-Mills theory with SU(2) gauge groupQ In this model 
the electromagnetic U(l) gauge symmetry is identified with the abelian U(l) subgroup of 
the SU{2). The other components of the SU(2) gauge field play the role of charged fields 
whose non-zero expectation values break the U(l) symmetry leading to a phase transition 
in the dual field theory. 

In this paper following [8] we will consider a five dimensional SU(2) Einstein- Yang-Mills 
theory with a negative cosmological constant whose action is given by^ 



S = / d b xy/=g 



-OR- A) - F a » u 



where F^ u is the field strength of the SU{2) gauge field. In our notation the negative 
cosmological constant is given by —12/L 2 with L — 1. The equations of motion coming from 



^ee also 11 for p-wave holographic superconductors in the context of five dimensional gauged super- 
gravity. 

2 In our notation we have set the five dimensional Newton constant to one, K5 = 1. 
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the above action are 



Rp, - \g,uR - Qg»„ = T^, -L^ t ( ^F^) + q r bc A\F c ^ = o, (1.2) 

^ V 9 



where 



T, v = Fl p F a v " - -g^F; v F a (1.3) 
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with F« v = d^A a v - d v Al + qf abc A b ^Al. 

The model has only one free parameter which in our notation is given by the gauge 
coupling q. Whether or not the back reactions of the gauge field on the metric are important 
is controlled by q. When q is large the effects of the gauge field on the geometry are negligible, 
while for finite q the gauge field back reactions are important and affect the geometry. 



The equations of motion (1.2) support an AdS solitonic solution with zero gauge field 



and the metric which is given by 

ds 2 = -^~-dr 2 + r 2 (-dt 2 + dx 2 + dy 2 )+r 2 g(r)d X 2 , g = \- T \. (1.4) 
r l g[r) r 4 

Essentially this solution can be obtained from a five dimensional AdS black hole solution by 
making use of two Wick rotations. Physically, this solution corresponds to a five dimensional 
solution compactified on a circle with anti-periodically boundary condition for the fermions 
along the compact direction. Since the geometry ends at r = r where g tt is non-zero, the 
background provides the gravity description of a three dimensional field theory with a mass 
gapQ 

It is worth mentioning that unlike the AdS black hole that finiteness of the gauge potential 
at the horizon prevents to have a non-zero A t at the horizon, in the present case one could 
still have the above solution with a constant non-zero gauge potential A t = fi. 



The theory given by the action (1.1) admits another analytic solution with non-zero gauge 
field. The corresponding solution is, indeed, an AdS Reissner-Nordstrom black hole which 
carries the charge of the £7(1) abelian subgroup of the SU{2) gauge group. The solution is 
given by 

ds 2 = — - gdt 2 + r 2 (dx 2 + dy 2 + dz 2 ), A = p ( 1 - a 3 dt, (1.5) 



where g = r 2 — 1+ r p 2 ^ + and a 3 is the generator of U(l) subgroup 



Here we have 



normalized the coordinates such that the horizon is located at r = 1. In this notation the 
Hawking temperature of the black hole is T = ^^Z 3 ' • 

The organization of this paper is as follows. In the next section we study p-wave holo- 
graphic insulator /superconductor in the probe limit where we will also compute the AC 
conductivity of the model. In section three we first consider fully back reacted solutions and 
then we study the phase structure of the theory. In section four, for completeness, we will 
study zero temperature limit of p-wave holographic superconductor in four dimensions. The 
last section is devoted to conclusions. 



3 We would like to thank referee for his/her comment on this point 
4 We denote the generators of SU{2) by a % for i = 1,2, 3. 
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2 Insulator/Superconductor transition at probe limit 



In this section we will consider dynamics of the SU(2) gauge field on the background(1.4). 
In general turning on a gauge field, the background metric (1.4) gets corrections due to back 
reactions of the gauge field back. Nevertheless at first order we will consider the case where 
the back reactions are negligible so that the gauge field may be treated as a probe. This can 
be done with the assumption that q is sufficiently large (9 10 . 



2.1 Phase transition 



To study the insulator/superconductor phase transition, following [9], we will consider the 
following ansatz for the gauge fielcj^] 



.4 



(r) a 3 dt + *0(r) a x dx. 



{2.1) 



Since we are in the probe limit, the back reactions of this gauge field on the metric (1.4) are 



negligible. Note that in this ansatz the t-component of the gauge field represents the U(l) 
gauge field, while the second term plays the role of the charged field whose condensation 
breaks the U(l) gauge symmetry. 



Plugging this ansatz into the equations of motion (1.2) one arrives at 



3 g' 
- + - 

r g 



2 



r 4g 



0. 



g_ 
g 



0/ + —-0 

r 4g 



0. 



(2.2) 



The aim is to solve the above equations, though in general it is difficult to do it analytically 
and therefore we will utilize the numerical method. To proceed we note that the above 
equations are invariant under the rescaling r — > r r, — > r o and ip — > r ip, so that r may 
be dropped from the equations and leaves us to work with dimensionless quantities. 
Near the boundary where r — > oo the solutions of the equations behave as 



li 



P 



0i 



(2.3) 



From dictionary of the AdS / CFT correspondence 12 13 , It is known that /i and p correspond 
to the chemical potential and charge density in the boundary theory, while ipo and ipi may 
be identified as a source and the expectation value of the dual operator, respectively. Since 
we are interested in the case where the dual operator is not sourced, we set ipo = 0. Also up 
to a normalization one has (O) ~ 

On the other hand imposing the finiteness condition near the tip where r — > 1, leads to 
the following series expansions for the ip and <p fields (see also [5]) 



a + a>i(l 



A> + ft(i-~)+ft(i-~) 2 



(2.4) 



5 See also 14 15 
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Figure 1: Here we plot the behaviors of (O) and charge density p for various values of 
chemical potential p. Here p c = 2.26. 



With these boundary conditions one may solve the equations (2.2) to find the expectation 
value of the dual operator as a function of the chemical potential. In fact given the initial 
values for the gauge field near the tip, one may find the boundary values of the gauge field 
by making use of the shooting method. Indeed doing so, one finds that the solution is 
unstable to develop a hair for the chemical potential bigger than a critical value, p > p c . 
More precisely using the Mathematica we have numerically solved the equations leading to 
the plot depicted in figure [TJ In our numerical method we have chosen <p at the tip as the 
shooting parameter and fixed ip at the tip by making use of the asymptotic condition x/jq = 0. 

From the figure [I] (right) we observe that the expectation value of the dual operator is 
non-zero for p > p c = 2.26. Having non-zero expectation value spontaneously breaks the 
£7(1) symmetry leading to a phase transition. To find the nature of the phase transition, it 
is illustrative to find the behavior of the charge density in terms of the chemical potential. 
Indeed since the charge density may be identified with the derivative of free energy, p = 
OF '/dp, its behavior around the critical chemical potential can indicate the nature of the 
phase transition. Indeed, as it is evident from figure [I] (right), it turns out that the phase 
transition is second order. In this sense the behavior of the system is similar to that for the 
s-wave one [H]. 



Since the geometry (1.4) is dual to a three dimensional theory with a mass gap, it is 
natural to treat the phase of the system before condensation as an insulator phase. Therefore 
the obtained phase transition may be interpreted as an insulator /superconductor phase 
transition. To see whether the new phase behaves as a superconductor it is useful to study 
the response of the theory to an external magnetic field. In other words one may compute 
the AC conductivity of the system. 

2.2 AC conductivity 

To study the conductivity of the theory we will consider an extra magnetic field along y 
direction. From gravity description point of view this can be done by turning on a non-zero 
gauge field in y direction. To do so, one may consider the following ansatz for the gauge 
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field along y direction 

A y = A{r)e iu)t a*. (2.5) 

Since we are working in the probe limit, we will consider the back reactions of the new 
component of the gauge field neither on the metric nor on the other components of the 
gauge field. As the result it is sufficient to consider the equation of motion for A y in the 



the equations (2.2). 



background generated by the metric (1.4) together with the solutions of ip and given by 



The corresponding equation of motion for A y is given by 

A" + ( - + A' + (u 2 -i, 2 )A = 0. (2.6) 
For large r where we approach the boundary the behavior of the gauge field is 

A = A + ^. (2.7) 



With this notation the conductivity in y direction is given by (see for example |10|) 

—iAi 

a w = -j-. (2.8) 

Using the numerical solution we had found for ip and <fi in the previous subsection we can find 
the behavior of the AC conductivity in terms of the energy co. Actually we expect that the 
imaginary part of the conductivity to have a non trivial behavior at u — 0. More precisely 
we would expect to get a delta function support for the real part of the conductivity at 
u) = 0. Indeed this is what we get from our numerical solution. 

In the figure [2] we have plotted the behavior of the conductivity in terms of the energy 
when the condensation is non- zero. At u> — > we find a pole showing that we have an infinite 
conductivity as expected for the superconductor phase. Note that the pole structure shown 
in figure [2] is repeated periodically. This is due to an infinite tower of vector modes. 



So far we have seen that the model based on the action (1.1) exhibits a p-wave insula- 



tor/superconductor phase transition as we vary the chemical potential. On the other hand 



as we have already mentioned in the introduction the equations of motion (1.2) admit an- 



other analytic solution; the RN black hole solution. Therefore one may expect that the AdS 



solitonic solution (1.4) could decay to an AdS black hole leading to another phase transition 
in the model. This is indeed the confinement /deconfinement first order phase transition 
considered in j6). In our context from three dimensional theory point of view this corre- 
sponds to an insulator /conductor phase transition. It is then natural to look for another 
conductor/superconductor phase transition in the model. 

We note, however, that to explore the whole picture of the phase structure of the model 
one has to consider the back reactions of the gauge field on the metric as well. In other 
words we will have to go beyond the probe limit. This is, in fact, what we will do in the 
next sectionB 

6 The complete phase diagrams for s-wave holographic insulator/superconductor system has been studied 
[7j where the authors have considered the back reactions of the gauge and scalar fields. 
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Figure 2: Here we plot the behavior of imaginary part of the conductivity , a yy , for various 
values of u when the expectation value of the dual operator is non-zero. Note that such a 
pattern repeated periodically. 



As a final remark, before going to the back reactions analysis, it is an instructive exercise 
to study the insulator/superconductor phase transition in the presence of a constant DC 
current (see for example (l6)). To do so we will consider the following ansatz for the gauge 
field 

A = (p(r)a 3 dt + ijj{r)a l dx + A(r)a 3 dy. (2.9) 

Note that although we are still working in the probe limit where the back reactions of the 
gauge field on the metric are negligible, all components of the gauge fields are treating in the 
same footing. This means that one should take in to account the effects of A y component 
on the other components of the gauge field. Indeed it turns out that the equation of motion 
for ip is modified as follows 

<+(* + ^V + £^ = 0, (2.10) 



while for <p it remains unchanged. Finally the last equation is given by (2.6) with u — 0. 

It is then the aim to solve these equations to find the behavior of the system. Actually 
we find that the model is again unstable against developing hair leading to the symmetry 
breaking when the chemical potential is bigger than a critical value. We note, however, that 
in the present case the critical value of the chemical increases as we increase the background 
DC current (see figure [3]). In other words having a non-zero background field makes it harder 
to have insulator/superconductor phase transition. 



3 Beyond probe limit 

So far we have been considering the model in the probe limit where the back reactions of the 
gauge field were negligible. The aim of this section is to study the effects of the gauge field on 
the background metric. We, note, that back reacted solutions of the five dimensional gravity 
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Figure 3: The behavior of the expectation value of the dual operator in terms of the chemical 
potential in the present of a constant background field Aq. The most left curve is for the 
case when A = 0, the others appear as we increase the field. Therefore as we turn on the 
background field, the condensation is harder to be formed. 



( 1.1 ) have been studied in [17] (see also [18]) where the authors have numerically constructed 
asymptotically AdS charged black holes of the model. It was also shown that, for sufficiently 
low temperature, the model develops vector hair. Generically this corresponds to a second 
order phase transition in the dual field theory. The aim of this section is to further study 



the back reacted solutions of the theory given by the action (1.1). 



3.1 AdS soliton 

In this subsection we will study the back reactions of the gauge field on the AdS solitonic 
solution considered in the previous section. To proceed, we start with the following ansatz 
for the metric and gauge field 

(If 

ds 2 = -^y + r 2 {-f{r)dt 2 + h{r)dx 2 + dy 2 ) + g{r)e~ x{r) drf , 
A = (j)(r)a*dt + i/)(r)a 1 dx. (3.1) 

Using this ansatz the Einstein equations of motion read 

/ / \r 2/ 2h 

h" ti /3 /' ti 

I + h 

h h\r 2/ 2h 



x" x [ — ~~ ~ — ~~~ I — — ir*9 — t*)" = — ■ — i — ■ — h 

\r 2g 2 J gr 2 r 2 h r 2 f r A fgh 

fh! (fh)' (I g' , . , . .. 

+ J 7T L [- + --x)- — = J1 r ! r^-, (3.2) 



fh fh \r g 
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r A fgh 
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2 y x'\ 1, 2 4x// V , <P' 2 , 3g 2 ^ 2 2 
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while from the equations of motion of gauge field one finds 



+ + ^ + 



r 2/ 2h g 2 J r 2 fg 

Now the aim is to solve the above equations though in general it requires numerical analysis. 
It is evident that the AdS solitonic solution considered in the previous section is, indeed, a 
solution of the above equation for tp = = 0. 

Near the boundary where r — » oo the behavior of the gauge field is as follows 

= /i-^, ^ = ^o + ^. (3.4) 

Of course we set 0>o = to avoid having source for the dual operator in the field theory. 
The main goal of this subsection is to find the behavior of ip\ as a function of the chemical 
potential fi for given electric charge q. For this purpose we utilize the numerical method 
using Mathematica. Essentially the procedure is very similar to what we have done in the 
previous section, though here we have to deal with more equations. 

To proceed we assume that at the tip, r = r , the function g goes to zero, i.e. g(r ) = 0. 
On the other hand one may set tq = 1 by making use of an evident scaling symmetry of the 
equations 

r ->• Ar, (t, x, y, 77) ->• A _1 (t, x, y, rj), g ->• X 2 g, (0, ip) X(i/j, ip). (3.5) 

Therefore by considering a series solution near the tip we are left with five parameters given 
/(l), h(l), x(l) and 0(1), -0(1). We note, however, that since the equations are symmetric 
under a constant shift in x, we may set x — at tip. On the other hand due to the scaling 
symmetries 

f^a 2 f, <p^a<p } (3.6) 

and 

h ->• b 2 h } V -> H, (3-7) 

one may fix f(l) and h(l) by setting, for example, / = h — 1 at the boundary. Altogether 
we are left with two parameters and q. Therefore for a given q we have a one parameter 
family of solutions parametrized by the value of ip at the tip. 

To summarize we can solve the equations by treating 0(1) as the shooting parameter for 
a fixed q and given Then we can scan the moduli space of parameters by changing q. 

Doing so we find the following picture. 

For generic q the behaviors of ip and p as the functions of the chemical potential are 
indeed the same as those we had in the probe limit (see figures [T]) . More precisely there is a 
critical chemical potential above which the system becomes unstable to develop vector hair 
leading to a second order phase transition in the dual field theory. 

We, note, however that in the present case the critical value of the chemical potential 
depends on the value of q. Indeed as we increase q the value of the chemical potential, where 
the phase transition occurs, decreases. An interesting phenomena we encounter is that there 



S 



is a critical value for q (q m 1.94) below which the system exhibits qualitatively different 
behaviors. This phenomena is also seen in the AdS charged black hole which we will discuss 
in the next subsection. Actually the situation is very similar to the s-wave holographic 
insulator/ superconductor system studied in (Tj. We will back to this point latter 



3.2 AdS charged black hole 



The model given by the action (1.1 ) admits AdS charged black holes when the back reactions 
of the gauge field are taken into account. These solutions have been numerically studied 
in 17 where the authors have also shown that the solutions are unstable to develop a vector 



hair for sufficiently low temperature. In this subsection, for completeness of our study, we 



will rederive the results of 17 . Of course our ansatz are slightly different from that in 17 



though our final results are similar them, up to a normalization. 



Motivated by 24 we consider the following ansatz 



ds 2 = 



dr 2 
g(r 



+ r 2 (h(r)dx 2 + dy 2 + dz 2 ) - g(r)e-^ r) dt 2 , 



(3.8) 



Plugging the above ansatz into the equations of motion (|1.2|) for thr metric one finds 

6 e 
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(3.9) 
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For finite q where the back reactions of the gauge field are important the model admits AdS 
charge black holes. In particular when we set ip = the equations admit an analytic solution 



which is indeed the Reissner- Nordstrom AdS black hole given by (1.5). 



For generic values of the gauge field, beside the RN black hole, there is no other well known 
analytic solution, though asymptotically AdS charged black holes have been numerically 
constructed in [17] (see also (l8)). Using the same procedure as that in the previous sections 
we may solve the equations numerically. To proceed we note that for large r as we approach 
the boundary the behavior of the gauge field is as follows 



P 

V ~ ^2- 



4> 



^° + ^2" • 



(3.10) 



Again we will set ^0 = 0. On the other hand near the horizon where g vanishes we assume 
<p is zero too. The horizon may be set to be located at r = 1. By making use of the 
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Figure 4: The behavior of the expectation value in terms of temperature for q = 3.5 where 
we get T c = 0.226. 



time reparamerization we can set \ to zero at the horizon. On the other hand utilizing the 
following symmetries 

r^Xr, g^X 2 g, (0, V) -> A(0, i/j), (3.11) 

and 

h->£h, (3.12) 

one may fix the values of g and h at the horizon by setting , for example, g = r 2 and h = 1 
at the boundary where we recove the AdS solution. 

Taking into account that the value of <fi at the horizon is treated as the shooting parameter 
we are left with two parameters given by q and ip at the horizon. Therefore for a fixed q we 
get one parameter family of solutions parametrized by the value of ip at the horizon. 

For generic q one finds that the solution is unstable to develop a vector hair for sufficiently 
low temperature where the U(l) gauge symmetry is also broken. This corresponds to a second 
order conductor /superconductor phase transition from field theory point of view (see figure 

3 (ED- 

We note that the critical temperature depends on the value of q. Indeed as we increase 
q the critical temperature increases too. It is worth mentioning that there is a critical value 
for q where the system shows qualitatively different behaviors. In the normalization of our 



numerical computations it happens at q ~ 1.94. This is exactly the behavior obtained in 17 
where the authors have found that there is a critical q where the transition becomes first 
order. This has to be compared with behavior we have found in the previous section. In the 
following section we will discuss about this point. 



3.3 Phase structure 

In this subsection we would like to explore the phase structure of the model using the 
results we have found so far. To proceed we note that the equations of motion support two 
distinctive solutions; AdS soliton and AdS charged black hole. In each case the solution 
becomes unstable to develop a vector hair as we change the parameters of the model. In 
particular we have seen that when we change the chemical potential there is a critical point 
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Figure 5: The qualitative phase structure of the model. As we change q the positions of 
points A and B are changed and therefore the shape of different regions get modified. There 
is a critical point where A and B are on top of each other. For sufficiently small q there 
is a possibility to have a first order phase transition where a superconductor becomes and 
insulator as we decrease temperature as shown in the left picture. The figure "a" is the 
phase structure in the probe limit, while the others are with back reactions. 



above which the AdS soliton develops a vector hair. From field theory point of view this 
phenomena corresponds to the insulator /superconductor second order phase transition. On 
the other hand when we have a charged black hole the solution is unstable to generate a 
vector hair for sufficiently low temperature which is indeed the holographic realization of the 
second order conductor/superconductor phase transition. 

On top of these we note that using the Euclidean solitonic solution one may associate 
a temperature to the solution which could be any value for a given chemical potential. 
Therefore as we draw the phase structure in the T — /x space the curve which separates 
insulator phase from the superconductor phase should be a line parallel to the T axes. On 
the other hand when the period of the Euclidean time of the AdS soliton becomes the same 
as that in the AdS charged black hole the free energy of the system shows that the favored 
solution is the AdS charged black hole and thus we have a first order phase transition from 
AdS soliton to AdS charged black hole [6]. So altogether we get four different phases as 
shown in figure [5}a. 

It should be mentioned that in order to find the above picture for the phase structure of 
the model the results we have gathered from probe limit considerations were sufficient. We 
note, however, that in order to fully explore different phases of the model in more detail one 
needs to go beyond the probe limit. In particular taking into account the back reactions we 
observe two new features in the phase diagram of the model. 

The first observation is that for large chemical potentail as we decrease temperature the 
favored phase is soliton superconductor. In other words the phase diagram gets modifed as 
figure [5}b. In particular we cannot have a phase describing hairy charged black hole at zero 
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Figure 6: The behavior of ha in terms q is shown by the dashed line, while the solid line is 
for fiB- Theses two lines meet each other at q « 1.94 where \xb ~ Ha ~ 1-87. 



temperature which could have been the case if the phase diagram had been given by figure 
[5]-a. We will back to this point in the next section. It is important to note that the figure 
[5] should be treated as a qualitative description and the interfaces between different phases 
have been ploted schematically. To explore the precise structure of the phase diagram higher 
numerical precision is needed, specially for small q. We hope to back to this point in our 
future work. 

On the other hand there are two special points, named by A and B in figure [5] which in 
order to understand their physical significant, it requires to study the system beyond the 
probe limit. 

Actually our numerical computations show that the positions of these two points labeled 
by ha and hb change as we are changeing q. More precisely, when we decrease q, ha and 
Hb increase as well and eventually they meat each other as we approach the critical point 
q = 1.94 where one gets Hb ~ Ha ~ 1-87 (see figure [6]) . 

As the result when we change the parameter q, different regions in the phase diagram of 
the model get modified and in particular we reach a situation where two points A and B are 
on top of each other as shown in figure [5]-c. 

This is indeed the critical point we have found in our numerical calculations in the 
previous subsections where the model shows qualitatively different behavior. This is also the 
situation which has been observed in [17] ( and also in (7j for s-wave) where it was shown 
that the phase transition becomes first order at the critical point. 

As we further decrease q we observe that while both ha and hb increase, the point A 
passes though the point B (see figure IHJ-d) . In this case we encounter a new phase transition. 
Actually as we decrease temperature there is a range of \i bewteen which the superconductor 
becomes an insulator via a first order phase transition. This phase transition has been also 
seen in the s-wave consideration in W\. Going further it seems that the phase where we 
have soliton superconductor becomes smaller and smaller and eventually disappears from 
the phase diagram, though due to the uncertainty of our numerical results, we have not been 
able to explore the situation exactly. In particular for low temperaute (for small enough q, 
i.e. q ~ 0.86) the numerical solution develops a singularity and one has to study T — > 
limit of hairy charged black hole more carefully. This is indeed what we will do in the next 
section. 
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4 Zero temperature limit 



So far we have been considering the phase transition between AdS soliton and AdS charged 
black hole with non-zero temperature. It is then natural to pose the equation what happens 
when we send the temperature of the AdS charged black hole to zero'iQ In general for a 
charged black hole sending temperature to zero we will end up with an extremal black hole 
whose near horizon geometry develops an AdS2 throat with non-zero entropy. Therefore if 
we would like to study holographic superconductors at zero temperature the extremal black 
hole cannot provide the gravity dual descriptions. In fact to get an eligible background one 
must have a geometry with zero size horizon ensuring that the ground state is a single state 
(entropy is zero). 

Indeed zero temperature limit of s-wave holographic superconductors in three dimensions 



has been studied in [19] (see also 20 23]) where it was shown that the corresponding hairy 
solution develops an AdS± geometry at the core of the space time. More precisely they have 
found numerical solutions which interpolate between hairy AdS± and AdS^ geometries. The 
Zero temperature limit of the holographic p-wave superconductors in three dimensions has 



also been studied in 24 



To complete our discussions of p-wave holographic superconductors in four dimensions, 
in this section we will consider zero temperature limit of the four dimensional p-wave holo- 



graphic superconductor. Actually to do this one needs to solve the equations of motion (3.9) 
with a particular boundary condition ensuring that the resultant geometry would have zero 
size horizon. To proceed we consider the following behaviors for the parameters of the ansatz 



03.9D in r ->■ 0+ limit. 

0~0 o (r), if) ~ Vo -ipi(r), X~ Xo-Xi(r), 9 ~ r 2 + #i(r), h ~ h + hx(r), (4.1) 

with the assumption that <po,ipi,Xij 9i an d hi go to zero sufficiently fast in the limit of 
r — > + . Plugging these behaviors into the corresponding equations of motion, at leading 
order, one finds 

= 0o^, X = Xo e r , g = r — — e r , 



* = *>( 1 - -4^ e ^)> h - h "{ 1 + ^^)- < 42 > 

where a = qipo/ho. On the other hand for large r one impose the following asymptotic 
conditions for the gauge field components 

^ = ^ - ^> ^ = ^ + ~k- ( 4 - 3 ) 

Now the aim is to solve numerically the equations of motion with the above boundary 
conditions. To do so, using the symmetries of the model we can set <p = h = 1 and Xo — 0. 

7 Note that the insulator/superconductor phase transition we have considered in the previous section was 
also occurred at zero temperature. We note, however, that this has not to be compared with the present 
case. This is because in the previous case to get the solitonic solution we have put anti-periodic boundary 
condition for the fermions along the compact direction which is not thr case here. 
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Figure 7: The behavior of <f),g/r 2 ,ip,x an d h as functions of r. The function is plotted for 
q = 0.93, 1, 1.12, 1.22, 1.32, 1.41 which have been shown by thick, dashed, dotted, red, blue 
and green curves respectively. 



On the other hand we would like to have a solution which is asymptotically AdS§ without 
hair. Therefore we are interested in the case where if) vanishes asymptotically. This can 
be done by imposing a relation between a and q. For example for 3 < q < 7 one finds 
a = 0.8 + 0.49g. 

It is straightforward to solve the equations numerically using Mathematica. The final 
results are shown in figure [7j As we see the solutions approach AdS^ geometries at IR with 
non-zero hair while they are asymptotically AdS*, without hair. Therefore from field theory 
point of view in IR we get an emergent CFT broken by the expectation value of the operator 
that is dual to if). 

It is worth noting that as we approach q — > v3/2 limit we lose the accuracy of our 
numerical solutions. This is the limit where we recover the extremal limit of the charged 
AdS black hole. Therefore effectively we will have to solve the equations on the AdS% 
geometry and y/3/2 is, indeed, the Breitenlohner-Freedman bound for mass in AdS2 space 
time. 

Having found hairy solutions it is natural to compute the conductivity for the model. To 
study the conductivity we consider small fluctuations for the metric and A y component of 
the gauge field as follows 

gty = 6/(r)e"^, A\ = 6^(r)e"^. (4.4) 

where e is a small number controls the effects of the perturbations. As far as the equations of 
motion for gauge field and metric are concerned, since we are working in the probe limit, the 
effects of the perturbations on the equations of motion of the other components are negligible 



and thus one can still use the equations (3.9). On the other hand for the perturbations, at 
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first order in e, the equations of motion are given by 

q 2 ijj 2 



I 1 -" • 

r g 

q 2 ijj 2 
r 2 gh 2 



x' ti 

— H 

2 /i 



0. 



e x cu 2 



g 2 r 2 gh 2 
r 2 f - 2rf + r 2 



+ 



J" + 6' 



1 /' x' h' 

- + — + — H 

r f 2 h 



0. 



(4.5) 



By making use of the equation of motion for 



(p +(p [- + 



g_ 
g 



2 



77 



(3.9), one arrives at 

q 2 i/j 2 



e x oj 2 



r2 



g 2 



r 2 gK' 



- e" 



9 



0. 



(4.6) 



To find the conductivity one needs to find the asymptotic behavior of gauge field (p. In fact 
for large r where we find <p = (po + ^ the conductivity is given by 



a 



-vpx 

LUlfiO 



(4.7) 



On the other hand at r — > 0, using the results of the last section and imposing in coming 
boundary condition at r — > we find 



(p(r -)■ 0) = ip 



6 r 



(4.8) 



where a; = \fu 2 — a 2 . 



Taking into account that the solutions we have found for the parameters of the anstaz 



(3.8) are real and using the equation (4.6) together with its complex conjugate we can show 
that the following expression is a constant of motion in the sense that it is r independent 

Im (rhge- x/2 (p*(p') . (4.9) 

Therefore its value at the horizon is the same as that at the boundary. Evaluating the above 
expression at the boundary and near r = one finds 2Im(^o^i) = ^ |<^o | 2 Re(d>) . Plugging 
this result into (4.7) we arrive at 
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Re (a) 



Re(cu) \(p 
2u \(p \ 



(4.10) 



that is zero for uj < a leading to a gap. Indeed, this is exactly (up to a factor of 1/2) the 
same expression which found in p24] for three dimensional model. 



5 Conclusions 

In this paper we have studied holographic insulator/superconductor by making use of a five 
dimensional Einstein- Yang- Mills theory with an SU(2) gauge field. The equations of motion 
supports both AdS soliton as well as AdS charge black hole solutions. Both solutions may 
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become unstable to forming a vector hair as we change the parameters of the solutions, 
including temperature and chemical potential. 

Therefore altogether the model exhibits four different phases. From gravity point of view 
they correspond to AdS soliton, hairy AdS soliton, AdS charged black hole and hairy AdS 
charged black hole and from field theory point of view they are insulator, superconductor, 
conductor and another superconductor, respectively. As we see there are two phases where 
the system is in the superconductor phase. Actually these two phases can be distinguished 
by the behavior of the conductivity [7]; while in the solitonic case the real part of the 
conductivity has a series of delta functions, in the charged black hole case we just have a 
gap at low frequency. 

Using the probe limit we have been able to study qualitatively the general features of 
the phase structure of the system, though to explore the phase diagram in more detail it 
was curtail to take into account the back reactions. Doing so, we have found that the model 
has rather a rich phase structure. In particular when q < 1.94 there is a range of \i between 
which the system can show a first order phase transition from superconductor to insulator as 
we decrease temperature. Indeed this is a new phenomena in which has been first observed 
in the context of s-wave insulator/superconductor transition in u\. It would be interesting 
to understand this phase transition better. 

We have also studied zero-temperature limit of four dimensional p-wave holographic 
superconductor. We have numerically solved the corresponding equations of motion have 
shown that equations of motion admit smooth solutions which interpolate between hairy 
AdS$ solution at IR and AdS?, solutions without hair at UV. Having found a hairy AdS$ 
solution at IR shows that there is an emergent CFT at IR such that the conformal symmetry 
is broken by the expectation value of the dual operator. Evaluating the conductivity for the 
solutions one finds that the real part of the conductivity vanishes for u < a leading to a gap 
in the model. 



Note: 

All results of the present paper are based on several Mathematica codes we have written for 
each part. We should admit that to develop our Mathematica codes the one prepared by C. 



P. Herzog which is available at his homepage ( http : / / wwwphy. princeton.edu / ^cpherzog /[ ) 
was illustrative. 
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